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Дается обоснование предлагаемых изменений в конструкции аварийного клапана химического реактора высокого давления. 
Численными методами теории упругости проведен анализ напряженно-деформированного состояния сопрягаемых деталей 
клапана. Параметрическими исследованиями получен вариант сопряжения элементов клапана, обеспечивающий герметизацию 
химического реактора. 


Введение 

Во время работы химического реактора рабочее 
давление газа в нем достигает 200 МПа. Аварийный 
клапан должен герметично закупоривать реактор 
до тех пор, пока давление в нем не превышает оп- 
ределенного предела, и сбрасывать излишки давле- 
ния в атмосферу, если рабочее давление превышает 
допустимое. Естественно, что аварийный клапан 

116 


должен иметь конструктивные особенности, кото- 
рые исключают утечку газа при штатном рабочем 
давлении. В действующей конструкции аварийно- 
го клапана между двумя сопрягаемыми металли- 
ческими поверхностями вставлялась серебряная 
проволока (кольцо), которая должна была служить 
герметиком (уплотнением) и предотвращать утечку 
газа через зазор между сопрягаемыми поверхностя- 
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ми. Недостатком такой конструкции является на- 
личие дополнительной сопрягаемой детали в узле, 
которая по регламенту подлежит замене при каж- 
дом ремонте или разборе клапана. 


1 



Рис. 1. Схема аварийного сбросного клапана: 1) штуцер ; 
2) нажимной фланец; 3) конусное седло; 4) корпус; 
5) канал; сплошной стрелкой показано направление 
движения газа в рабочем режиме, штриховой при 
сбросе давления; 6) плунжер, стрелкой показаны 
направления его движения 

Наряду с действующей предлагается более 
простая конструкция клапана (рис. 1), в которой 
сопрягаемые металлические поверхности имеют 
вид усеченного конуса. Преимущество такого тех- 
нического решения заключается в устранении до- 
полнительной детали в узле сопряжения и упроще- 
ние процедуры сборки клапана при его ремонте. 
Если угол конусности седла клапана и основания 
клапана (корпуса) совпадают, напряжения в зоне 
контакта равномерно распределены по всей пло- 
щади усеченного конуса. Если же угол расточки 
седла клапана и его основания отличаются на один 
или несколько градусов, то напряжения в зоне кон- 
такта возрастают, а сама зона контакта при несов- 
падении углов конусности уменьшается. 

Величина зоны контакта и параметры напря- 
женно-деформированного состояния (НДС) седла 
и основания клапана рассчитываются исходя из со- 
отношений теории упругости. 

Основные соотношения теории упругости. 

Метод решения 

В осесимметричной задаче теории упругости 
кручение отсутствует, компонента перемещений ѵ 
вдоль координаты Ѳ равна нулю, компоненты ми» 
не зависят от Ѳ (рис. 2) [1]. 

Деформации вдоль осей г, Ѳ, і, а также дефор- 
мация сдвига в плоскости нравны: 

_дУ _Ц_ ш _дЦ д\Ѵ 
8 "~ дг’ 8ѳѳ ~ г’ 8г2 ~ дг дг ' 

Ненулевые компоненты тензора напряжений 
связаны с компонентами тензора деформаций сле- 
дующими соотношениями: 
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или в компактной форме: 

{ст}=[0]{8}. 



Рис. 2. Оси и перемещения вдоль осей цилиндрической сис- 
темы координат 


Осесимметричная задача теории упругости в 
данной работе решается методом конечных эле- 
ментов (МКЭ) [2]. Компоненты перемещений II, 
]Ѵ аппроксимируются внутри треугольного конеч- 
ного элемента линейной функцией. Соотношения 
между деформациями и перемещениями в матрич- 
ной форме имеют вид: 

{8} = [В]{11}. 


[5] - матрица градиентов - имеет вид: 
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где 


= г і 2 к - V, ; ъ і = 2 і - 2 к ; с, = Г к - г к ; 


^■ = 2 ^■(а,-+6,г + с і г) ; А = ± 


г к 


Коэффициенты а р Ь р с р а к , Ь к , с к функции формы 
Щ, И к получаются круговой перестановкой индек- 
сов і, і, к. 

После того, как определены матрицы [Б] и [В], 
нетрудно определить матрицу жесткости конечно- 
го элемента для осесимметричной задачи теории 
упругости в виде: 

[К е } = \[В] т [0][В]сІѴ, (1) 
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где (ІѴ=2лгсІА элементарный объем тора, получен- 
ного вращением сечения сІА вокруг оси. 

Если в формуле (1) матрицы [П] и [5] содержат толь- 
ко постоянные величины, то они могут быть вынесены 
из под знака интеграла. Однако матрица [5] содержит 
коэффициенты, являющиеся функциями координат. 
Заменим переменные величины г и іш средними зна- 
чениями 7 ѵй каждом конечном элементе, тогда матри- 
ца жесткости для конечных элементов в осесимметрич- 
ной задаче теории упругости примет вид [2]: 

[К е ] = ][В] [0][В]-2 л.гА, 

V 

где [В\ - матрица коэффициентов, в которой вмес- 
то переменных величин г, і используются их сред- 
ние значения по элементу. 

Матричное уравнение для ансамбля элементов 
имеет вид: 

[К]-{Щ = {Р} 

где [.ЭД -глобальная матрица жесткости, которая соби- 
рается из матриц жесткости элементов отдельных эле- 
ментов [К ] , {7} - глобальный вектор-столбец нагрузки. 

Вектор-столбец нагрузки {/) для ансамбля эле- 
ментов формируется следующим образом. Пусть на 
тело вращения воздействует поверхностная нагруз- 
ка р, и р-. Тогда вектор узловых нагрузок 

(2 к,+к ] ) Рі : 

(2К,+К і )р : 

(Я І +2Я ] )р г 

(Я, +2 Я,)Р; ’ 

0 
0 

где - длина стороны между узлами / и у. 

Если рассматривается горизонтальная поверх- 
ность, то и тогда на наиболее удаленный от оси 
вращения узел будет приходиться большая часть 
нагрузки, чем на узел, расположенный ближе к оси 
вращения. Если рассматривается вертикальная по- 
верхность, то К=Ц и компоненты нагрузки поровну 
распределены между узлами конечного элемента. 

Для получения разрешающей системы уравне- 
ний МКЭ проводится процесс ансамблирования 
конечных элементов по всей расчетной области и 
получается глобальная матрица жесткости и гло- 
бальный вектор нагрузки для всей области. 

Система линейных алгебраических уравнений, 
получающаяся в результате применения процеду- 
ры МКЭ, симметрична, содержит 938 уравнений, 
имеет ленточную структуру (ширина ленты равна 
40), решается методом Гаусса [3]. 

Реализация метода пошагового нагружения 

Расчетная область, которая включает в себя 
контактирующие поверхности аварийного клапа- 
на, изображена на рис. 3. 


На участках контура по нормали к поверхности 
действуют давления: о^І^ЮОО МПа; а^^ОО МПа; 
ст ггІге“-330 МПа. Касательные напряжения на этих 
площадках принимаются равными нулю. 

Перемещения на участках СО и ОН равны нулю: 
ѵ|се=0, Ч се =0, ѵ| ея =0, Ч ея =0. 


г 



Рис. 3. Конфигурация расчетной области. Многоугольник 
АВСйЕЕ — седло клапана, многоугольник ССНОР — 
основание клапана, линия 00 — ось вращения 

На рис. 3 угол конусности седла клапана равен 
28°, угол конусности основания клапана равен 30°, 
и угол рассогласования составляет 2°. 
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Рис. 4. Сетка конечных элементов 

Седло клапана и основание клапана могут 
иметь различные углы конусности, поэтому сетка 
конечных элементов наносится независимо на сед- 
ло клапана и на основание клапана (рис. 4). 

Особый интерес представляет изменение сетки 
в зоне контакта седла клапана и основания клапа- 
на в процессе нагружения. Пусть в начальный мо- 
мент нагружения (рис. 5) площадь контакта седла и 
основания клапана равна длине стороны одного 
элемента, умноженной на 2л г. Тогда в начальный 
момент приложения нагрузки узлы 86 и 85, 105 и 
104 связаны (совпадают), а узлы 123 и 124 находят- 
ся на некотором расстоянии друг от друга. По мере 
возрастания сжимающей нагрузки седло и основа- 
ние клапана деформируются и узлы 123 и 124 пос- 
тепенно приближаются и сливаются в один. После 
этого седло и клапан имеют общую площадку, рас- 
положенную между узлами 86-105-124 с одной сто- 
роны на седле и узлами 85-104-123 с другой сторо- 
ны на основании клапана. 



118 






Технические науки 


106 124 



Рис. 5. Сетка конечных элементов в зоне контакта 


Задача расчета клапана решается методом пошаго- 
вого нагружения [4]. На каждом шаге процесса рас- 
считывается НДС седла и основания. Затем контро- 
лируется, насколько близко подходит узел на седле 
клапана (узел 124) к соответствующему узлу на осно- 
вании клапана (узел 123). В процессе решения на каж- 
дом шаге происходит вычисление узлов деформиро- 
ванной сетки конечных элементов. Если расстояние 
между соседними узлами седла и основания клапана 
меньше наперед заданной малой величины, узлы свя- 
зываются в один. Процесс нагружения может продол- 
жаться дальше, однако площадь контакта на следую- 
щих шагах нагружения увеличивается по сравнению с 
начальной. Таким образом, площадь контакта седла и 
основания клапана по мере развития процесса нагру- 


жения будет изменяться, что и учитывается в алгорит- 
ме последовательного нагружения. 

Процедура связывания узлов состоит в следую- 
щем. Перемещения в связанных узлах равны, следо- 
вательно, строки и столбцы в глобальной матрице 
жесткости, соответствующие перемещениям связы- 
ваемых узлов, не являются независимыми. Один из 
связываемых узлов принимается за основной, дру- 
гой за вспомогательный, строки и столбцы вспомо- 
гательного узла складываются со строками и столб- 
цами основного узла. Затем строки и столбцы вспо- 
могательного узла преобразуются следующим обра- 
зом: диагональный член умножается на ІО 6 , а компо- 
ненты глобального вектора нагрузки, соответствую- 
щие вспомогательному узлу, приравниваются нулю. 

После решения системы алгебраических урав- 
нений получается, что перемещения во вспомога- 
тельном узле равны нулю. Приравнивая их переме- 
щениям в основном узле, получаем реальную кар- 
тину перемещений в упругом теле. 

Расчет напряженно-деформированного состояния 
аварийного клапана при различных углах 
конусности седла 

Сначала рассчитывается НДС узла седло-основа- 
ние клапана в случае, когда угол конусности седла и 
основания клапана совпадают и равны 30°. Поверх- 
ности и линии уровня интенсивности напряжений и 
деформаций в случае одинаковых углов конусности 
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Рис. 6. Поверхности интенсивности напряжений <т, (а) и интенсивности деформаций е ; (с), линии уровня интенсивности напря- 
жений (Ь) и линии уровня интенсивности деформаций (б). Углы конусности седла и основания клапана совпадают (30°) 
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Рис. 7. Поверхности интенсивности напряжений а ■ (а) и интенсивности деформаций е ; (с), линии уровня интенсивности нап- 
ряжений (Ь) и линии уровня интенсивности деформаций (сі). Углы конусности седла (28°) и основания клапана (30°) 
не совпадают 


приведены на рис. 6. Параметры НДС клапана в слу- 
чае, когда углы конусности седла (28°) и основания 
клапана (30°) не совпадают, приведены на рис. 7. 

Следует отметить, что параметры НДС, рассчи- 
танные для изменения углов конусности седла в 
пределах от 27° до 29°, отличаются незначительно. 

На рис. 7 видна локализация параметров НДС на 
небольшом участке контактной поверхности при не- 
совпадении углов конусности седла и основания 
клапана. Угол конусности седла клапана 28° обеспе- 
чивает на поверхности контакта седла и основания 
клапана уровень напряжений, сравнимый с уровнем 
напряжений на поверхности контакта седла и плун- 
жера. Так как при этом уровне напряжений на конта- 
ктной поверхности "седло-плунжер" обеспечивается 
герметизация, то и на поверхности "седло-основание 
клапана" обеспечивается герметичность соединения. 

Таким образом, в конструкции аварийного кла- 
пана, в котором угол конусности седла клапана и 
основания отличаются на 1 ...3°, обеспечивается 
герметизация на поверхности "седло-основание 
клапана" за счет локализации сжимающих напря- 


жений, сравнимых с напряжениями на контактной 
поверхности плунжер-седло клапана. 

Выводы 

1. В работе приведен пример численного решения 
осесимметричной контактной задачи теории уп- 
ругости. Особенность постановки задачи и ее ре- 
ализации заключается в том, что область контак- 
та заранее неизвестна, а ее размеры определяют- 
ся и уточняются в ходе решения. Разработанный 
алгоритм, основанный на использовании МКЭ в 
сочетании с процедурой пошагового приложения 
нагрузки, позволяет получить оценки параметров 
НДС в зоне контакта на каждом шаге нагружения 
и в окончательном рабочем состоянии. 

2. Анализ НДС сопрягаемых элементов клапана в 
зоне контакта показывает, что изменение гео- 
метрии сопрягаемых элементов приводит к тако- 
му уровню напряжений в зоне контакта, что это 
сопряжение становится герметичным. В связи с 
этим предлагаемые упрощения в конструкции 
клапана являются вполне обоснованными. 
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